Collapse dynamics of a ^^^Yb-^^^Yb Bose-Einstein condensate 
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In this paper, we present a theoretical study of a two-component Bose-Einstein condensate com- 
posed of Ytterbium (Yb) isotopes in a three dimensional anisotropic harmonic potential. The 
condensate consists of a mixture of ^^^Yb atoms which have a negative s-wave scattering length 
and ^^^Yb atoms having a positive s-wave scattering length. We study the ground state as well as 
dynamic properties of this two-component condensate. Due to the attractive interactions between 
^^"^Yb atoms, the condensate of ^^^Yb undergo a collapse when the particle number exceed a critical 
value. The critical number and the collapse dynamics are modified due to the presence of ^^*Yb 
atoms. We use coupled two-component Gross-Pitaevskii equations to study the collapse dynam- 
^ , ics. The theoretical results obtained are in reasonable agreement with the experimental results of 

^ ■ Fukuhara et al. [PRA79, 021601(R) (2009)]. 

PACS numbers: 03.75.Hh, 03.75.Kk, 03.75.Mn 

^ , I. INTRODUCTION 

The first experimental observation of Bose-Einstein condensate (BEC)i"— in bose atom vapors have initiated an 
^ ' exciting field of research, both theoretically and experimentally. One of the most interesting developments in this 
field is the formation of multi-component condensates. Multi-component BECs have been observed experimentally 
qh" by Myatt et alA and Hall et al^ for two different hyperfine spin sates of ^''Rb, by Modugno et al^ for different 
^ ; , atoms (^^K and *^Rb), and Papp et al^ for different isotopes of the same atom (^^Rb and ^^Rb). A rich variety 
d • of various interesting effects exhibited by these two-component BECs have inspired a number of theoretical studies 
^ ] covering various aspects of a these systems^^— . The common feature of the bose systems in these experiments is 
that the intra-componcnt and the inter-component boson-boson interactions are all repulsive. It rises the curiosity 
' O I about the properties of a multi-component condensate in which one kind of atoms have repulsive interactions while 
another kind of atoms have attractive interactions. Recently, Fukuhara et ali^ observed BEC of spin-zero Yb isotopes 
by implementing an all-optical cooling protocol. The bose-bose mixture in these experiments contain ^^^Yb atoms 
having a positive s-wave scattering length and ^^^Yb atoms having a negative s-wave scattering lengt h^^d^ . The s- 
wave scattering length between ^^"'Yb and ^^^Yb is also positivei^. Such a two-component condensate can be expected 
, to show dynamical properties far more complex than a one-component condensate of attractively interacting bosons, 
' which becomes unstable when the number of atoms exceed a critical value^ d^'^^ . 
ly-j [ For a two-component BEC with attractive interactions between bosons in one component and repulsive interactions 
QQ I in the second component, two most basic questions are of that of its stability and that of the collapse dynamics. In this 
' paper, we theoretically study the ground state as well as dynamic properties of such a two-component condensate in 
\^ . a three dimensional (3D) anisotropic harmonic potential. We specifically consider the ^^^Yb-^^^Yb bose-bose mixture 
(•~^ ' in the anisotropic harmonic confining potential as in the experiment of Fukuhara et ali^. The paper is arranged as 
] follows. In Sec. II, we describe the theoretical model for the study of a two-component BEC. In Sec. Ill, we discuss 
^-H • the ground state profile of a stable two-component BEC. In Sec. IV, we present the collapse dynamics of the system 
^ I and compare it with the experimental resultai^. The conclusions arc given in Sec. V. 
, 

H , II. TWO-COMPONENT BEC: THEORETICAL MODEL 

The ground state and dynamic properties of a two component Bose-Einstein condensate (BEC)is well described by 
a set of coupled Gross-Pitacveskii equationai^"— (GPE) given by 

(r t) 

ih ' — +v,ir)+gumr,t)\^ +g,,\ij,ir,t)\^)ij,{r,t), (1) 

where i = 1, 2 are indices for the two components (1 for ^^^Yb and 2 for ^^^Yb), j ~ 3 — r = {x, y, is the spatial 
coordinate vector, Vi{r) = [1 / 2)mi{uy^x'^ -^ux^y^ -l-w^z^) is the trapping potential, gn = ATrh^au/mi is the intra-species 
interaction and gij = 2Trh'^ aij / rriij the inter-species interaction strength between atoms in the condensed state, an is 
intra-species and aij is inter-species s-wave scattering length, rriij = minij / [rrii + mj) is the reduced mass in which 
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FIG. 1: Ground state profiles of ^^''Yb (top panel) and ^^®Yb (bottom panel) condensates in a 3d anisotropic harmonic potential 
with A^i = 6 X 10'' and N2 = 150. 

mi and rrij are atomic masses. The normahzation condition for each component is j \%j:i{r)\'^dr = Ni , where Ni is 
number of atoms in each component. 

We non-dimensionalize Eq. ([T]) through a set of linear transformations: t = LOxt, f = r//, ^/'i(r) = ^/•^"-0(r). 
After dropping the wiggles on the symbols, we obtain 



i^^t^^ = (-iv^ + v,{v) + \umv,t)\' + A.,|V^,(r,t)|^)V.(r,t) , (2) 

where 



l — \ l , Aii — ; , Aij — ; , 

mUJx l ' 

^,(r) = i(x2 + ^V+7V), ^^.'^^ ^='±. (3) 

Since the masses of the two isotopes are nearly equal, we have taken mi — m2 = m. In order to find a stationary 
solution of Eq. (I2|), we do a separation of variables ipiijc^t) = tpiir) x cxp[—i{ij,i/{hidx))t], where is the chemical 
potential of the ith component. Starting from Eq. we obtain 

(-iv2 + i;,(r) + XMr)\^ + X^,\tP,{r)\^)Mr) - -^Mr) • (4) 



III. GROUND STATE PROFILES OF A TWO-COMPONENT BEG OF YB ATOMS 



In this section, we discuss the ground state properties of a two-component BEC of Yb isotopes by numerically 
solving the coupled GPE (Eq.H]). The ground state solution of the GPE is found by the imaginary time propagation 
method. In this method, the time dependent GPE is evolved in imaginary time starting from an initial guess using 
a finite difference Crank-Nicholson (FDCN) scheme^!. In imaginary time propagation we have taken the space step 
as 5x = 5y = 5z ^ 0.1 and the time step as 5t = 0.00005. We have used a set of parameters corresponding to the 
i74Yb-i76Vb system in the experiment^: m = 2.8734238 x IQ'"^^ Kg, an = 5.55 x 10"^ m, 022 = -1.28 x IQ-^m, 
ai2 = 021 = 2.88 x 10"^ TO, yx{=0Jx/'2.iT) = 4:5 Hz, {= ujy/2'K) = 200 Hz, v^{^ojj2-n) = 2,00 Hz. 

Due to the attractive interaction between ^^^Yb atoms, the condensate of ^^^Yb undergo a collapse if the particle 
number exceeds a critical value N2ir^'^. For a single component ^^^Yb this value is given by ~ 0.5L/|a|; where 
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FIG. 2: Time evolution of number of '^^*Yb(filled circles), '^^^Yb (open squares, open circles, filled squares). The corresponding 
experimental results are given by (open triangles for ^'^^Yb and filled triangles for ^^®Yb). The values of three body recombina- 
tion terms are; K3 = 4.2x 10"^^cm''s"^(for all points)), Ki = 3.0 x 10"^** 0771*^3"^ (open squares), K3 = 3.0 x 10"^'' cm'^s"^ (open 
circles), Ki = 3.0 x lO^^^cjn^s"^ (filled squares). The time is in units of seconds. 



L = ^yh/mui and u! = (wj^Wj^wz)^ . For the parameters given above, A^2c ~ 250. But, the critical number is modified 
due to the presence of ^'''*Yb atoms with a positive scattering length. The physical origin of this modification is in 
the effective potential felt by the bosons in the attractive component. The interaction with the repulsive component 
changes the effective potential from V2{r) — |A22||'02(r)P to V2{'r) — |A22||V'2(r)|^ + X2i\ipi{r)\'^ . This reduces the 
effect of the attractive interaction and leads to a reduction of the critical number. Alternatively, one may say that 
the mean-field contribution from the repulsive component leads to a flattening of the effective potential felt by the 
attractive component. If A^i = 6 x 10'', then for the given parameters, A*2c is calculated to be ~ 220 by our numerical 
simulation. In order to prepare a stable condensate, we have taken A^i = 6 x 10^ and N2 = 150. The ground sate 
profile of the two-component BEC is presented in Fig. 1. We observe that the ^''^Yb atoms are at the center of trap 
and are surrounded by a large cloud of ^'''^Yb atoms. 



IV. COLLAPSE DYNAMICS OF A TWO-COMPONENT BEC OF YB ATOMS 

In this section, we study the collapse dynamics of a two component BEC composed of '^^^Yb and '^^^Yb atoms using 
the coupled time-dependent GPE's. Due to the negative s-wave scattering length, the -'''^Yb condensate becomes 
unstable if the number of atoms becomes greater than a critical value N2c- This condensate collapses by emitting 
atoms out of it. Due to high density of atoms in the attractive condensate, the loss of atoms from the condensate 
occurs through three-body collisions. To model this collapse we add a imaginary three body quintic loss tern^"— to 
the RHS of GPE Eq. ^ given by 

if.V'. = -(l/12Ki^^|V.|V,, (5) 

where is the three-body loss coefficient for each component. We have neglected the two-body dipolar loss term as 
they make negligible contribution in this casei^^"— . We have also left out the loss terms for '^^^Yb — ^^"^ Yb — ^'^^ Yb, 
i74Yb -176 Yb -17*^ Yb, i^^Yb _i76 Yb -174 Yb, and ™Yh -"^ Yb -"^ Yb collisions, since the losses due to these 
arc comparatively smalli^. So, Eq. ([2]) becomes 

»^^^ = (-^V2 + 7;,(r) + A,,|^,(r,t)|2 + A,,|V^,(r,t)|2-ze,|V',(r,i)|')V'.(r,t), (6) 

where C» = il/l2)N^K'2l-^Lj-\ 

We time evolve the coupled time-dependent GP equations Eq. (|6]) using the finite difference Crank-Nicholson 
(FDCN) scheme with a known initial condition. In real time propagation, we have taken the space step as 
5x = Sy = 6z = 0.1 and the time step as St = 0.005. The time evolution of the number of i™Yb and ^^^Yb is 
shown in Fig. 2 along with the experimental data of Fukuhara et alt^. Considering the complex dynamics of the 
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FIG. 3: Ground state profile of ^'^''Yb at different times. Tiie four figures correspond to t{seconds) — 0.0 (top left panel), 0.2 
(top right) 0.4 (bottom left), and 2 (bottom right). The values of three body recombination terms are: A'l = 4.2 X 10"^'' cm s , 
Kl = 3.0 X IQ-"^^ cm^s-^. Here, A/'i = 6 x IQ-* and N2=2x 10*. 



the two-component system with mixed interactions, the theoretical results may be said to be in reasonable agreement 
with the experimental results. We observe that there is a significant loss of ^^^Yb atoms. We also see that the decay 
of ^^^Yb is very rapid. It is due to the collapse of the ^^^Yb condensate. The number of ^^^Yb atoms does show a 
very small decrease, which is not visible on the scale of this figure. To understand the details of the decay process, 
we study the condensate profiles of each component at different times. The results are shown in the Fig. 3 for ^^^Yb 
and in Fig. 4 for ^^^Yb. At t = 0, the ^^^Yb are at the center of the trap (top left panel of Fig. 3) surrounded by 
^^■^Yb (top left panel of Fig. 4). Since the number of atoms in the attractive component is higher than the critical 
number for stability, the system is unstable. When the system evolves in time, the attractive component explodes as 
is evident from the spreading of this component with time, in real space, as shown in Fig. 3. The spiky structures in 
these figures represent the inhomogcneities produced due to the on-going explosion process. As mentioned earlier, the 
explosion also leads to a spread of the ground state profile. Due to the coupling between the attractive and repulsive 
components, the condensate of ^''"'Yb is also redistributed in real space during the time evolution, as shown in Fig. 
4. The numbers of remaining atoms in each condensate component during the time evolution, for a longer period of 
time, is shown in Fig. 5. We note that the best agreement with the experimental results (Fig. 2) is obtained for the 
measured values of the K\ and K^. The disagreement at later times is likely to be originating from the neglect of 
atomic loss due to collisions involving ^^^Yb and ^^''Yb. During the initial stages of the time evolution, the bosons 
distribution gets heavily mixed due to the explosion and due to the coupling between the two components. Then, at 
later times, the inter-component collisions is likely to affect the atom loss. We are unable to include these loss terms 
since their values are not known at present. 



V. CONCLUSIONS 



In this paper, we presented a study of some static and dynamic properties of a two-component bose condensate 
consisting of repulsively interacting ^^^Yb atoms and attractively interacting ^^^Yb atoms in an anisotropic harmonic 
confinement. In the stable state, the ground state has ^^^Yb atoms in the center of the trap surrounded by ^^^Yb 
atoms. When the number of atoms in the attractive component exceed a critical value, the the system undergo a 
collapse. We analyzed the time evolution of this collapse process for the specific system parameters of the ^^^Yb- 
^^"'Yb system studied in the experiment of Fukuhara et ai. The details of the collapse dynamics are found to be 
in reasonable agreement with experimental results. The critical number for stability of the attractive condensate is 
reduced by it's interaction with the repulsive condensate. 
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FIG. 4: Ground state profile of ^^*Yb at different times. Tiie four figures correspond to t{seconds) — 0.0 (top left panel), 0.2 (top 
right) 0.4 (bottom left), and 2 (bottom right). The values of the three body recombination terms are: = 4.2 x W~^^cm^s~^ , 
Kl = 3.0 X IQ-'^^cm^s-^. Here, A/'i = 6 x lO^* and N2=2x 10*. 
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FIG. 5: Time evolution of number of ^^*Yb(top curve), ^'^^Yb(solid line, dashed, dash-dot lines). The values of three body 
recombination terms are: Kl = 4.2 x 10"^'' cm'^s"^ (for all lines)), A'| = 3.0 x 10"^* cm'^s"^ (solid line), Kl = 3.0 x 10"^'' cm^s"^ 
(dashed line), Kl = 3.0 x 10"^'' cm^'s"^ (dash-dot line). The time is in units of seconds. 
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